Decoherence mechanisms of 209 Bi donor electron spins in isotopically pure 28 Si 
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Bismuth ( 209 Bi) is the deepest Group V donor in silicon and possesses the most extreme char- 
acteristics such as a 9/2 nuclear spin and a 1.5 GHz hyperfine coupling. These lead to several 
potential advantages for a Si:Bi donor electron spin qubit compared to the more common phospho- 
rus donor. Previous studies on Si:Bi have been performed using natural silicon where linewidths 
and electron spin coherence times are limited by the presence of 29 Si impurities. Here we describe 
electron spin resonance (ESR) and electron nuclear double resonance (ENDOR) studies on 209 Bi in 
isotopically pure 28 Si. ESR and ENDOR linewidths, transition probabilities and coherence times 
are understood in terms of the spin Hamiltonian parameters showing a dependence on field and 
mi of the 209 Bi nuclear spin. We explore various decoherence mechanisms applicable to the donor 
electron spin, measuring coherence times up to 700 ms at 1.7 K at X-band, comparable with 28 Si:P. 
The coherence times we measure follow closely the calculated field-sensitivity of the transition fre- 
quency, providing a strong motivation to explore 'clock' transitions where coherence lifetimes could 
be further enhanced. 



Amongst the first proposals for quantum information 
processing (QIP) in solid state devices was that by Kane 
in 1998, using phosphorus (P) dopants in silicon, the first 
of the group V donors pQ. The choice of the P-dopant 
in particular can be attributed to its use in classical sili- 
con technology as donors, in addition to the simplicity of 
the spin system which consists of an electron and nuclear 
spin 1/2. Si:P has been extensively studied as a potential 
qubit 2 , 3 J and is more than ever one of the leading can- 
didates in the quest for a solid-state quantum computer 
[U[5]. It begs the question, however, of whether other 
donors of the same group, i.e. arsenic (As), antimony 
(Sb) and bismuth (Bi) , could have the same if not better 
properties than P. As the deepest group V donor, Bi has 
the largest nuclear spin (9/2) and the largest hyperfine 
coupling (1.4754 GHz) [B], and has received increasing 
attention over the past two years [71-fTT]. A single donor 
has a 20 dimensional Hilbert space, allowing either more 
storage space or a more robust encoding of information. 
The large hyperfine coupling enables shorter nuclear ma- 
nipulation times and a zero field splitting of 7.3 GHz, 
making Si:Bi useful as memory for hybrid superconduct- 
ing circuits [T2J [13]. At low magnetic fields, the spin 
Hamiltonian provides so-called 'clock' transitions for the 
electron spin, where the transition frequency between two 
states is insensitive to first order in magnetic field fluc- 
tuations of the environment [M14T6] . On the other hand, 
high vapour pressure and low solubility of the donor in 
silicon complicates the doping process during growth [T7] , 



while the high atomic weight (209) increases defect den- 
sity during ion implantation [llj . 

Previous works [THlOj on bismuth were realized using 
doped natural silicon crystals where the presence of other 
silicon isotopes, in particular 29 Si with a nuclear spin of 
1/2, resulted in spectral broadening and coherence times 
limited to about 1 ms. We investigate here isotopically 
pure 28 Si crystals [IS], where the Bi was introduced dur- 
ing crystal growth using the method developed earlier 
for natural Si [IT]. The residual concentration of 29 Si is 
50 ppm, allowing us to explore the intrinsic properties of 
the bismuth donor electron spin in silicon. Pulsed elec- 
tron spin resonance (ESR) experiments were performed 
using a X-band (9.74 GHz) Bruker Elexsys E680 and 
E580 spectrometers. Unless otherwise mentioned, the 
concentration of activated Bi is 2 x 10 15 cm~ 3 . 

The bismuth donor coupled nuclear/electron spin sys- 
tem is described in frequency units by an isotropic spin 
Hamiltonian: 

Hu = B { lE S z ®t+ lN t®I z )+AS.T (1) 

where the two first terms correspond to the electronic 
(S) and nuclear (I) Zeeman interactions with an external 
field Bo and the last term corresponds to the hyperfine 
coupling A (Fig. [I]). In addition to the normal ESR- 
allowed transitions [Arris = ±1; Am/ = 0] (high field no- 
tation), eight 'forbidden' transitions [Ams — ±1, Am/ = 
+2] can be found at the lower field side of the ESR 
lines for -7/2 < rrtj < +7/2 (Figure^). The separa- 
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FIG. 1: Eigenenergies of Si:Bi in frequency units. Twenty 
eigenstates arise (for Bo > 0) from the electron spin S — 1/2 
and the nuclear spin / = 9/2. Eigenstates are coloured ac- 
cording to their mi value in limit of large Bo, and allowed 
transitions at X-band (9.7 GHz) are highlighted. Inset shows 
the corresponding X-band electron spin echo-detected spec- 
trum. 



tion between allowed and forbidden peaks increases from 
-30 /zT for uij = +7/2 to -270 /zT for to/ = -7/2. Sim- 
ulations using EasySpin [19] confirm the non-negligible 
intensity of these forbidden transitions, especially for low 
| to/ | values which are the most mixed in terms of ms 
and to/ . An electron spin echo (ESE) intensity ratio 
between the allowed/forbidden transitions of nearly 0.25 
was measured after optimization of the microwave power 
to each. Therefore, these forbidden transitions offer an 
additional way to manipulate more states of the 20 di- 
mensional Hilbert space at the speed of a typical ESR 
pulse (tens of ns). 

These |Ato/| = 2 transitions were not observed in ear- 
lier studies using nat Si:Bi because unresolved hypcrfine 
coupling to the surrounding 29 Si broadens the ESR lines 
to about 0.4 mT. In contrast, the continuous wave (CW)- 
ESR linewidths in 28 Si:Bi vary between 8.2 /zT and 22 /zT 
(Figure^). Their to/ dependence can be explained by 
a full width half maximum (FWHM) hyperfine distribu- 
tion between 40 and 90 kHz (depending on the sample), 
and a static B field inhomogeneity about 3 /zT (a typi- 
cal value for ESR). The hyperfine spread corresponds to 
about 30 ppm of the total value and arises most likely 
from lattice strains (the crystal was found to contain 
dislocations). Figure [5J: shows electron- nuclear double 
resonance (ENDOR) linewidths measured for about half 
of the 36 transitions. These can be explained by the 
same distribution in the hyperfine coupling A (the Bq 
inhomogeneity has negligible influence on the ENDOR 
linewidth) . 

The electron spin coherence time Ti a was measured us- 
ing a Hahn echo to eliminate the effect of static inhomo- 
geneities. Using this method, a of 20 ms can be mea- 
sured at 5 K, however this technique adds its own type of 
decoherence to the system known as instantaneous diffu- 
sion (ID). The microwave 7r-pulse in the tt/2 — n — echo 
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FIG. 2: a) Typical allowed (Am/ = 0) and forbidden tran- 
sitions (| Ami | = 2) observed by ESE-detected field sweep. 
The microwave power is optimized separately for each tran- 
sition. b,c) Simulated ESR and ENDOR linewidths in light 
grey are calculated assuming a Gaussian spread in A and Bo 
of about 60 kHz and 3 fiT respectively. Empty points (model) 
and filled points (experiment) correspond to the transitions 
at 9.74 GHz. ESR and ENDOR linewidths were measured at 
25 K and 10 K respectively. For the ENDOR linewidths, the 
colours are guides to the eye to facilitate comparison of the 
experimental and simulated values. At the lowest measured 
field, the ENDOR linewidths only provide an upper bound 
due to experimental limitations. 



sequence flips all the spins in the sample, creating a 
change in the local magnetic field experienced by each 
spin. This additional dephasing mechanism, purely ex- 
perimental, can be reduced by shortening the 7r-flip to a 
smaller tip-angle rotation, at the cost of less signal/noise. 
By extrapolation to a tip angle of zero, the intrinsic T2 e ,i n t 
can be found 1201: 



1 

?2c 



C2o,int 



+ C ■ (2ir leS ) 2 



(2) 



where 6 is the tip-angle of the second microwave pulse 
and C is the effective spin concentration (l/10 th of the 
actual donor concentration as only one out of the 10 
ESR lines is excited by a given pulse and participates 
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FIG. 3: Measurements of Ibe.int for the ten ESR transitions 
in Si:Bi. a) ID slopes resulting from the fit of T^e times ob- 
tained using 7r, 7r/2,7r/4 and tt/8 pulses as measured at 5 K 
for each mj (some examples are shown in inset). The slopes 
vary according to (df/dB) 2 , as expected from dipole-dipole 
interactions between neighbouring Bi electron spins, b) Ex- 
trapolated T2e,int from the fit in (a) in the limit where 6 — > 0. 
The model shown in blue takes into account various sources 
of fluctuating magnetic field. 



in ID). 7 c ff is the effective gyromagnetic ratio, equal to 
df/dB, the gradient of the transition frequency with re- 
spect to the magnetic field. It represents both the sen- 
sitivity of the spin to the environment and its magnetic 
dipole strength [21] . 

For each sample with a given Bi donor concentration, 
we measured T 2c as a function of tip-angle for each of 
the different to/ (inset in Fig. [3^,), and fit the results to 
Eq. |2]) with a single concentration C and 7 c fr for each 
to/ as fitting parameters. In all cases, the best-fit value 
of C matched that extracted from electrical resistivity 
measurements to within 10%. The fitted values for 7 c g- 
are shown in Fig. [3^, and are in excellent agreement with 
the theoretical df/dB for each transition, as calculated 
from the spin Hamiltonian in Eq. (1). In addition to the 
slope, the intercepts at zero tip-angle were also extracted 
from the fit providing the extrapolated T^int which are 
summarized in Figs. [3] and [4] as a function of to/, tem- 
perature and donor concentration. 

The different values for df/dB for the ten ESR lines of 
Si:Bi offer an ideal opportunity to explore the contribu- 



tions of decoherence mechanisms for donor electron spins 
in silicon, as studied previously in 28 Si:P [2]. T^e is lim- 
ited directly by T le above about 9 K, while below this 
temperature decoherence is driven by spin flips (T lc ) or 
energy- conserving flip-flops of neighbouring donor elec- 
tron spins. Under both of these mechanisms, the effect 
on decoherence of a central spin is the product of both 
the sensitivity df/dB of the central spin and the effec- 
tive fluctuating magnetic field produced by neighbouring 
spins. This effective field is related to the magnetic dipole 
strength of the neighbours, df/dB, however, since each 
neighbouring donor can be in any one of ten to/ states, 
the relevant df / dB value must be averaged over different 
to/ at the particular value of the applied magnetic field 
B . Therefore, the effect of (T le ) spin flips of neighbours 
can be simply modelled by an average dipole strength 
(df/dB). On the other hand, the effect of flip-flops be- 
tween neighbours is the average product of the flip-flop 
rate and the magnetic field strength. At X-band, calcu- 
lations show that this contribution has a similar trend to 
Tic of neighbours. In summary, the overall decoherence 
rate can be expressed as: 



1/T 2 , 



df_ 
dB 



P \dB 



(3) 



where a is an additional field fluctuation, possibly from 
the external magnet, and /? represents both neighbour- 
Tie and flip-flop mechanisms as they are indistinguishable 
here. This is in excellent agreement with the observed 
experimental to/ dependence of T 2 e,int, shown in Fig. [3]} 
at 5 K. 

The temperature dependence of T 2Ci i nt in Fig. [4^, arises 
from the combination of the (temperature independent) 
neighbour flip-flop process and the (strongly tempera- 
ture dependent) T\ e . The latter appears twice in the 
sum of rates which determine T 2c : directly as T\ c for 
spin flips of the central spin and as y/Ti e for spin flips 
of the neighbours [2]. T le is independent of concentra- 
tion for small donor concentrations [3T] . As temperature 
decreases, Ti e increases above seconds and neighbour flip- 
flops become the limiting factor for coherence. 

By fitting the temperature dependence of T 2c ac- 
cording to the mechanisms above, we extract a low- 
temperature limit of T 2c corresponding to the neighbour 
flip-flop mechanism, and plot this time in Fig. |4Jd against 
concentration for the three samples we have studied here. 
For comparison, 28 Si:P values from Ref [2] under similar 
experimental conditions are also shown. 

For concentrations above a few 10 14 cm -3 , the effect 
of flip-flops play a limiting role in T 2c , which follows a 
quadratic dependence on concentration 22 . This de- 
pendence can be understood by considering that concen- 
tration affects both the flip-flop rate, and the sensitivity 
of the central spin to a neighbour flip-flop. At lower con- 
centrations, the effect of residual 29 Si impurities becomes 
non-negligible and limits the coherence time to about 2 s. 
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FIG. 4: Measurements of extrapolated T2 e ,int as a function 
of temperature and donor concentration, a) Temperature de- 
pendence of T2e,int and Tie for different donor concentrations. 
As temperature is lowered, T2 e ,int is limited first by 2\ e di- 
rectly, then by spin flips (Tie) of neighbouring Bi electron 
spins, and finally by flip-flops of the neighbours. T lc was mea- 
sured using saturation or inversion recovery. Around 6 and 
8 K, stretched exponential decays were observed as spin flips 
(Tie) of neighbours becomes the dominant mechanism. In 
these cases, a better measure of the coherence is the "mem- 
ory time", shown as T2 C here, when the curve's amplitude 
reaches 1/e. b) Fitted decay times in the limit where T — > 
are shown from (a) for 28 Si:Bi, and from Ref [2] for 28 Si:P. The 
donor concentration dependence is expected to be quadratic 
for flip-flops until the effect of residual 29 Si nuclear spins be- 
comes dominant at very low donor concentration. 



The temperature dependence of T% e for Si:Bi was stud- 
ied by Belli et al. [10] and found to be driven by an Or- 
bach process above 25 K. Below this temperature, a first- 
order T -7 Raman process has been suggested in recent 
papers [SrfTU]. however we find that a T~ 9 Raman pro- 
cess from the spin-orbit interaction to be a better fit. 
The former process (T~ 7 ) is expected to be frequency- 
dependent (Tic oc / _2 ) while the latter (T~ 9 ) has no 
frequency dependence [33] ■ We measured T lo as a func- 
tion of temperature at 7 GHz and found the same values 
as at 9.74 GHz, following T" 9 down to 5 K. 

We have measured a T20 time of Bi donors in 28 Si of 
up to 700 ms at 1.7 K (using a refocussing pulse of an- 
gle 0.087T radians), comparable with that for P donors. 
Throughout this study of the various decoherence mech- 
anisms present in 28 Si:Bi we have observed that the elec- 
tron spin decoherence rates follows df / dB to at least first 
order. This provides a motivation to study T2 C at various 
'clock' transitions present in the range 5.2-7.3 GHz, at 
which df I dB drops to zero. At these points, we would ex- 
pect to have a greatly reduced sensitivity to donor or 
29 Si concentration, and approach the fundamental limit 
of Tic which exceeds one second below 6 K. 
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The difference in sensitivity to 29 Si spectral diffusion be- 
tween Bi and P can be accounted for, within experimental 
errors, by the difference in df / dB (a factor of 0.65). 

The average dipole coupling rate between two near- 
est neighbour Bi donor electron spins is about 10 Hz for 
a donor concentration of 10 15 cm -3 . In contrast, the 
ESR linewidth is due to inhomogeneity in both Bq and 
A corresponding to broadening of about 60 and 200 kHz, 
respectively. Both of these linewidths are much greater 
than the dipole coupling rate and would largely suppress 
any nearest neighbour flip-flops if they arose from in- 
homogeneity across a length-scale comparable to donor- 
donor separation. For this reasons we conclude that ob- 
served spread in hyperfine coupling strength A derives 
from long-distance variations during the crystal growth 
process, consistent with strain fields from dislocations. 
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